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ABSTRACT 

The equilibrium statistical mechanics of a d dimensional "oriented" manifold in an N + d dimensional 
random medium are analyzed in d = 4 — e dimensions. For N = 1, this problem describes an in- 
terface pinned by impurities. For d = 1, the model becomes identical to the directed polymer in a 
random medium. Here, we generalize the functional rcnormalization group method used previously to 
study the interface problem, and extract the behavior in the double limit e small and N large, finding 
non-analytic corrections in 1/N. For short-range disorder, the interface width scales as u; ~ Z>, with 
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. We also analyze the behavior for disorder with long- 
range correlations, as is appropriate for interfaces in random field systems, and study the crossover between 
the two regimes. 



PACS numbers: 75.60.Ch, 75.10.Nr, 64.60.Ak 



I. Introduction and Summary 

Oriented elastic manifolds embedded in spaces which contain random impurities that can pin the man- 
ifold occur in many physical systems. Both the dimension of the manifold, d, and the dimension of the 
space in which it is embedded, D = d + N , where N is the number of transverse dimensions, play important 
roles. [1,2] The case d = 1 corresponds to a directed polymer in a random potential [3, 4] which describes the 
interaction of a single flux line in a type-II superconductor with impurities. [5] Interfaces between two coex- 
isting phases in D dimensional systems are D — 1 dimensional oriented manifolds whose properties control 
much of the behavior of such systems in the presence of randomness. [6-8] It has also been argued that flux 
lattices in superconductors have an intermediate distance regime in which they behave like a d — N = 3 
oriented manifold. [9] 

From a theoretical point of view, it has become clear in the last few years that elastic manifolds in 
random media exhibit much of the interesting and subtle phenomena that characterize spin-glasses and 
other complicated disorder dominated phases. [4,10] In addition, the equilibrium statistical mechanics of 
directed polymers in a random potential can be mapped to the dynamics of an interface growing by random 
deposition - leading to insights into both problems. [11] In this paper we will analyze the equilibrium behavior 
of manifolds in random media for d just below the critical dimension of four, focusing on the limit of large 
N. 

By definition, an oriented manifold has no overhangs; therefore, it can be described entirely in terms of 
a set of N transverse coordinates {4> 1 }, which are functions of the d internal coordinates {x J } parameterizing 
the manifold. [12] The Hamiltonian describing such a manifold in the presence of a quenched random potential 
F(0(x),x)is 

H = Jd d x {iv0-V</> + V(0(x),x).}, (1.1) 

where 0(x) £ is a vector describing the transverse coordinates of the manifold with internal coordinate 
x £ » d . 

The random potential and/or thermal fluctuations will generally cause the manifold to roughen, resulting 
in divergent fluctuations in <fi. It is conventional to parameterize these by a roughness exponent (: 

(([^(x)-^(x')] 2 ) t )^|x-x'|^, (1.2) 

where the inner angular brackets with a T subscript denote a thermal average, while the outer brackets 
indicate a statistical average over the random potential. 

In the absence of the random potential, the interface will be flat for d > 2, corresponding to £ = 0, but 
be thermally rough with £ = £t = (2 — d)/2 for d < 2 (with logarithmic corrections in two dimensions). For 
d < 2, weak randomness is irrelevant when N > Nt — 2d/(2 — d), and an unpinned phase exists at high 
temperature. [13] For a strong random potential, low temperature or outside this regime of the N — d plane, 
the randomness always dominates and the system is controlled by a non- trivial zero-temperature fixed point. 

Much is known about the case of one internal dimension (d = 1), owing both to the mapping to interface 
growth[ll] and to the simplicity of numerical simulations. [14] When both d = 1 and N = 1, the roughness 
exponent is known to be exactly £ = 2/3. [3,15] For TV > 2, in addition to the low temperature randomness- 
dominated phase, a high-temperature phase emerges, in which the disorder is irrelevant and £ = 1/2. The 
value of £ in the low-temperature phase has been investigated numerically for TV > 2; it decreases with N 
and appears to approach 1/2 as N — > oo.[14] The possibility of a finite upper critical dimension, such that 
£ = 1/2 for N > N c , has been suggested by several authors [2, 16], but no evidence of this has appeared in 
rather extensive numerical simulations, and others have argued that no such upper critical dimension exists, 
but rather that £ decreases continuously to 1/2 as N — > oo.[4] 

For manifolds with d > 4, a perturbative analysis[8] shows that the interface remains flat (£ = 0). In 
this regime perturbation theory, or equivalently a simple RG treatment (briefly described in section II) is 
valid. For d < 2, as mentioned above, both high temperature pinned and low-temperature randomness- 
dominated phases exist for large enough N. Between two and four dimensions, only the disorder-dominated 
phase exists, characterized by a non-trivial £. The first attempts to analyze this phase involved perturbative 
methods about four dimensions - e.g. Parisi and Sourlas [17] - which yield incorrect results (in particular 
£ = (4 — d)/2) due to the existence of many extrema of the Hamiltonian Eq.(l.l). In order to carry out a 
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proper e expansion, a functional renormalization group (RG) is needed. This was introduced by Fisher [8] 
for the case N = 1, yielding the results £ = e/3 and £ = 0.2083e for random-field and random-bond disorder 
respectively to lowest order in e = 4 — d. The first result [7] is believed to be exact for 1 < d < 4, while the 
second is an O(e) calculation, which required a numerical solution of the RG fixed point equation. 

More recently, Mezard and Parisi[l] (MP) have performed an approximate "replica symmetry breaking" 
calculation on the model for general N, and argued that their results are exact in the limit as N — > oo. In 
MP's method, replica symmetry breaking is introduced as a variational ansatz. An interesting question to 
address is whether the replica symmetry breaking corresponds to some physical aspect of the problem (at 
least at large N) or is merely a feature of the restricted variational ansatz. This is a particularly interesting 
issue in light of our O(e) treatment, which does not involve replica symmetry breaking, but is a systematic 
pcrturbative RG calculation. It is hoped that a comparison of the two methods may provide insight into 
whether "replica symmetry breaking" has any well-defined meaning. 

In this paper, we generalize the 4 — e RG calculation of ref.[8] to arbitrary N. In the limit of large N, 
the fixed-point and stability equations become tractable, and we perform an expansion around this limit, 
working always only to first order in e. For short-range correlated disorder, we find 
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(1.3) 



while for a random potential with long-range correlations transverse to the manifold, 



(V(4>, y)VW, x'))c ~ |0 - 0'r <5 d (x - x'), 



(1.4) 



there is a continuously variable exponent which is independent of N 



CLK-e/(4 + 7), 



(1.5) 



in agreement with previous results for the particularly interesting case of interfaces in random-field systems 
which corresponds to N = 1 and 7 = — 1. The long-range fixed point becomes unstable when (lr < (sr- 
As in the N = 1 case, the long-range result is believed to be exact, while the short-range result is true only 
to 0(e). 

Before proceeding with the RG calculation, we briefly outline the remainder of the paper. In the next 
section, the model is described, and an RG procedure is developed to analyze the zero temperature fixed 
point in 4 — e dimensions. It is shown that a perturbative expansion, which incorrectly deals with the physics 
of the many mctastable states, breaks down. Section III analyzes the correct solution for the fixed point 
and roughness exponent in the large N limit, in which analytic results can be obtained. The stability of this 
fixed point is analyzed in section IV. The behavior for long-range correlated disorder (e.g. random fields) is 
studied in section V, as well as the stability of the associated fixed points. In section VI we summarize our 
conclusions and open questions and suggest several possible directions for future work. Appendices A and B 
contain various technical details, while Appendix C rederives the RG relations by an iterative minimization 
of the Hamiltonian and discusses the appearance of many minima. Finally, Appendix D analyzes possible 
multicritical fixed points. 
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II. Model and Zero- Temperature Renormalization Group 

The partition function in terms of the Hamiltonian Eq.(l.l) is 

Z{V} = J[dt/>] expt-if/T). (2.1) 
We take the random potential V to have a Gaussian distribution, with the two point correlation function 

(V(<f>, x)V (0', x')) = R{</> - <t>')5 (d) (x - x'), (2.2) 

with a short-distance scale implicitly included in the <5-function. To order e = 4 — d higher cumulants will 
be shown to be irrelevant. 

To organize a renormalization group (RG) treatment of the problem, we employ the replica method of 
averaging over the disorder. Note, however, that we do this only to organize the perturbation expansion. 
The partition function is now a random variable, and rather than follow the flow of its distribution function 
directly under the RG, one can follow the full set of moments. These are readily averaged over, yielding 
terms of the form 

W = j[d<t>] cxp(-ff p ), (2.3) 



with 



P =H^£ V0 Q • V0«- ^ RW-tP) - ^3 E S(<p a -^ a -^ -■], (2.4) 

J y a a/3 a/3j J 

where the terms with three or more replicas result from non-Gaussian correlations in the disorder; these 
are generated at higher order in e. We have not included terms with additional gradients, which will 
be irrelevant. The symmetry under simultaneous shifts of all the replica coordinates corresponds to the 
statistical symmetry of the disorder under spatial translation normal to the manifold. We have used p 
instead of the more conventional n for the number of replica indices to avoid confusion with the transverse 
dimension N of the manifold. 

The momentum shell RG approach we employ consists of integrating out high-momentum modes in a 
shell with A/6 < |p| < A. We take the limit in which the width of this shell is infinitesimal, i.e. b = e dl , 
which simplifies the formulae somewhat. To keep the cut-off fixed, momenta, coordinates, and fields are 
rescaled according to 

P = P% 

x = 6x', (2.5) 
0(x) = btytf). 

A simple first attempt at an RG analysis consists of expanding the function R(<j>) in a power series, and 
analyzing the results term by term. Since R is an even function, only even powers of (f) appear in such an 
expansion, i.e. 

a,f3 m 

Above four dimensions, the quartic and higher vertices are irrelevant and the theory flows to a simple 
Gaussian fixed point. Below four dimensions the quartic interaction becomes relevant, and one might hope 
to make a simple epsilon expansion by going to second order. A simple calculation shows that no perturbative 
fixed point exists, so that a strong-coupling analysis is necessary. 
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Physically, the behavior below four dimensions is dominated by the randomness, and should be described 
by a zero temperature fixed point. By allowing temperature to renormalize, it is possible to organize an 
analysis of this fixed point. The RG flows which arise from the scale changes are then 



dT 

Hi 

dR(4>) 



= (2-d-2C)T, 

sc 



dl 



dl 



sc 



sc 



(4-d-4C)i?(0)+C^ii(0), 

(6 - 2d - 205(0!, 4> 2 ) + Ci^udu + faidxWfa, 4> 2 ), 



The C,(j> 1 diR(4>) term in the second equation comes from the infinitesimal field rescaling of Eq.(2.5). 

The flow equation (Eq.(2.7)) for the temperature is actually exact, due to the "Galilean" invariance of 
the Hamiltonian Eq.(l.l). This can be seen by considering the change in the free energy corresponding to 
a uniform tilt of the manifold (or equivalcntly a change in boundary conditions). If the fields are shifted 
by a linear function of the coordinates cf> — > <f> + v l x l , the probability distribution of the new free energy is 
identical to that of the old one plus an additive constant 

AF = A(-TlnZ) = ^\v\ 2 L d . (2.8) 

This is an exact statement about the model, and must therefore be true at all stages of the RG; it requires 
that T (i.e. the coefficient of the stiffness term) only be renormalized by the scale changes. [18] We thus have 

f = -*r, (2-9) 

with = 2( + d — 2 determining the scaling of energies at the zero temperature fixed point. 

One can again attempt to proceed by expanding the function R(<p) in a power series. The condition 
fixing the quadratic term R 2 in Eq.(2.6) is then £ = (4 — d)/2 = e/2. However, since £ is of order e, all the 
higher terms in the expansion of R{<$>) also become relevant below four dimensions. It is therefore necessary 
to keep track of the entire series of {R m }, or the entire function R(4>). 

The one-loop RG equations are best derived in their functional form, through the use of Fourier trans- 
forms. As an example, consider the first order feedback of the R(<j>) term. As noted above, it will not 
renormalize itself, but it will contribute to the free energy. The first step in the calculation is to Fourier 
transform the interaction term, 



-Him = i E / e" 6 -^*)-^*)) (2.10) 

21 a j3 K 

where we use Greek and Latin letters for Fourier transforms perpendicular (<fi direction) and parallel (x 
direction) to the manifold, respectively, and f xK = f d d xd d K / '(2Tr) d . 

To perform the elimination part of the RG transformation, the fields are split into slowly and rapidly 
moving parts, 

0(x) = <Mx)+<Mx), (2-H) 

and a trace is performed over the "fast" fields </> > (x). In Eq.(2.10), the slow and fast terms separate into 
two exponential factors, and it is a simple matter to average over the fast modes. Dropping the label (<) 
on the slow fields, the traced term becomes 

-(Hint) =2^E/ e iK <^^-^^e-^ TG >^- s ^R{K), (2.12) 
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where the free two-point function is 

G>(x) = r e iqx M (2-13) 

where the > on the integral denotes integration over momenta (d d q/ (27r) d ) in the shell only. When evaluated 
at x = 0, this function has the well-defined limit 

G>(0) = (2ir)- d S d A d - 2 dl = A d dl, (2.14) 




where Sd is the surface area of a unit sphere in d dimensions. The second exponential in Eq.(2.12) can thus 
be expanded to yield 

'H int ) = 4^ / l Y, - 4f) - pd&IKp) } , (2.15) 



where the di act on the internal coordinates of the functions R{4>) (not on the spatial coordinates x) and 
repeated Latin indices are summed from 1 to N. The first term is of the appropriate form to feed back into 
R(4>), but it is reduced by a factor of T from the term in Eq.(2.4). Physically, this renormalization is due to 
the averaging of the potential by thermal fluctuations, and is thus negligible at the zero temperature fixed 
point of interest. The second term in Eq.(2.15) contributes to the renormalization of the free energy density 

f = (d - 6)f - ^d t d t R(0) + ■■■. (2.16) 

We now proceed with the analysis to second order in R(<p); one must multiply two terms and take the 
connected expectation value over </>>. Performing the average, one finds the expression 



E / exp{iK 1 -(^(x 1 )-^(x 1 ))+i K2 .(^(x 2 )-^(x 2 ))} 

(2.17) 



g T 4 

Oil, 01 xi/t 
"2, 02 X2,K 2 



x exp {-«?TG>(0)(1 - S aifSl ) - ^TG>(0)(1 - S a202 )} R( Ki )R(k 2 ) 
x exp • k 2 TG>(xi - x 2 ) [<5 Q1Q2 + 8 plfh - S aif>2 - S^ 2 } } 
— Disconnected Parts. 



The terms resulting from expanding the G>(0) parts are canceled by disconnected pieces, leaving only the 
expansion of the final term. Expansion of the final exponential gives terms proportional to 1/T 3 and lower 
powers of T. The 1/T 3 term is a three-replica contribution, 

^E / i ^(f(x)-f(x))a,]J(f(x')-/(x'))G > (x-x'), (2.18) 

since a single Kronecker delta function leaves three free replica indices. Eq.(2.18), however, contributes only 
at large momenta, since G> (x) contains an integral only over momenta in the shell. Because the generated 
3-replica term exists only at large momentum, it cannot feed back to generate a zero-momentum term until 
second order. The resulting contribution to the renormalization of R will then turn out to be higher order in 
e, since the 3-replica piece will be 0(e 2 ), anticipating that R will be O(e) at the fixed point. The situation 
is analogous to the neglect (to lowest order) of the (j> 6 term in momentum-shell RG for the conventional X(j) 
theory, in which a </> 6 interaction with large momentum is generated from two <j> terms, but does not feed 
back in a dangerous manner. Fig.[l] shows diagrammatically how the three-replica term is generated and 
feeds back into the two-replica piece. 
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The 1/T 2 parts generate both two- and three-replica terms. The three- replica term is down by a factor 
of T, and can be neglected at the zero temperature fixed point. Keeping track of the factors, one finds a 
contribution 



dR(cj>) 



dl 



o(R 2 ) 2T2 a ,fs 



/ G >( x - x ') 2 k9 3 fl(fW-/(x))9 ! 9 3 fi(fM-/(x')) 
a ^ ' X ' L (2.19) 

-29 l a J i?(</) a (x)-0 /3 (x))9 l a J i?(o) . 



Since we arc interested in the rcnormalization of the long-wavelength portion, the kernel K(yi) = G>(x) 2 
must be evaluated at zero momentum, i.e. 



/MP = 01= / fMxF= /' exp^(p + pO-x] = j> l_ , S;|2t L f// (2 2()) 



x.p.p' P'P" Jp P^ 



By rescaling R((f>) by a constant multiple, this factor can be removed. The full RG equation for R then 
becomes 



=(4 - d - 4C)i?(0) + C^dM) + 
+ 0(R 3 ). 



d l d j R{(j))d i d 1 R{(f)) - ^^^(0)^9^(0) 



(2.21) 



This equation can also be formally derived by expanding R{4>) in (fi > and keeping terms up to O(0>), 
although such a treatment does not properly treat the case of non-analytic R((f>). The flow equation Eq.(2.21) 
is exactly equivalent to the infinite series of RG equations obtained from all one loop diagrams in an ordinary 
diagrammatic approach. It was used previously in Ref.[2] and Ref.[19]. In Appendix C, we derive the RG 
equation schematically by directly minimizing H over the fast degrees of freedom > without the use of 
replicas or field theoretic techniques. 

At this point, one may try to directly analyze the flows and fixed points of Eq.(2.21): if a fixed point 
R <~ e is found, then the other terms in H p will not play a role to O(e). We will take this approach in the 
next section, but first it is instructive to investigate the manner in which the standard polynomial RG breaks 
down. A conventional RG approach is equivalent to expanding R{4>) in a power series 

and following the flow equations for the coefficients. The first two equations are 

--{e-2QR 2 , 



dR 2 



dk N + 8 (2 ' 23) 

__ £jR4 + __ i?4 . 

While the quadratic term could be fixed by requiring Q = e/2 (the naive perturbative result), R4 then flows off 
and is not stabilized at second order. Examination of the flow equation for R4 shows that it becomes infinite 
after a finite amount of renormalization. This is an artifact of the truncation to second order; nevertheless, 
i?4 will rapidly become O(l). 

We thus see that there is no fixed point of the perturbative RG for which R{4>) is analytic. It is the 
assumption of analyticity which leads to this conclusion, and we shall see that non-analytic fixed points of 
Eq.(2.21) do exist. The correct behavior for small 4> can be found by directly examining Eq.(2.21) with 
dR/dl = 0. This will be done in the next section. 

Before proceeding, it is useful to recall how a normal perturbative analysis leads to drastically wrong 
results. Various methods have been used to show that, in a field theoretic expansion, formally R 2 has no 
non-trivial renormalizations.[17] An apparent fixed point can then be found by setting ( = e/2 to all orders 
in e. This results in the so-called "dimensional reduction" result 9 = 2. In the supersymmetric formulation 
of Parisi and Sourlas which averages over all of the extrema of H with positive and negative weights, the 
only operator which appears corresponds to R 2 . Thus this analysis completely misses the flow of the other 
operators, such as R4, out of the regime in which a perturbative analysis might have been valid. As we 
discuss in the conclusion, it is likely that there is an exact upper bound for £ which is violated by the naive 
perturbative result. 
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III. Fixed Point Analysis in the Large TV Limit 

In this section, we analyze the behavior of the fixed points in 4 — e dimensions in the limit of large N. 
We look for solutions where the R is 0(l/N) by rescaling R — ► R/N. The flow equation (Eq.(2.21)) then 
takes the form 



9R(4>) = (4 _ d _ 4C)i?( 0) + (^dM) ' ! 



IdidjRWdidjRty) - didjRWdidjRQ) 



dl v " * r " 7 N . 

For rotationally invariant solutions, the ansatz = Q((j) 2 /2) yields 



(3.1) 



^ = (4 - d - 4C)Q + 2(yQ' + l -{Q'f - Q'Q' + 1 [2yQ'Q" + 2y 2 (Q") 2 - 2yQ"Qy , (3.2) 
where primes denote differentiation with respect to y, defined as 

V = ^ 2 /2, (3.3) 

and 

Qo = Q'(y = o). (3.4) 

Differentiating once and pulling out the 0(e) factor by defining 

Q' = -(^-d-2()-^—u, , 

+ 2 ' (3.5) 

and t=(4-d-2Ql, 

gives the final form of the fixed point equation: 

du 
~dt 

with 



= u + l3yu' + u'(u(0) -u)-fj, [3y(u') 2 + 2y 2 u'u" + yu"{u - u(0))] = 0, (3.6) 



f3= 2C/(4-d-2C), 
H= 2/(N + 2). 



(3.7) 



Any physical solution of this equation will have a finite value of u at the origin, so that by a choice of scale 
we can set 

it(0) = 1. (3.8) 

The behavior at large y will fix the value of (3 in a way which is somewhat analogous to more conventional 
eigenvalue problems. If we start with potential correlations which are short range, it is natural to expect the 
fixed points also to correspond to short range correlations. We thus look for fixed points with R((f>) decaying 
rapidly for large </>; i.e. u(y) decaying rapidly for large y. 

Naively in the large N limit the terms proportional to fi in Eq.(3.6) can be simply dropped, and the 
resulting problem has a much simpler form: 

u oo +pyu' oo +u' oo (l-u oo ) = 0. (3.9) 

We will see that this approximation is not valid globally but this equation determines the primary solution, 
called the outer solution in the language of boundarydayer theory. [20] This is in a sense the solution for 
^ = 00. The boundarydayer, or region in which this solution is not valid, occurs for large y. The solution 
in this tail region is conventionally denoted the inner solution (as it is valid inside the boundary layer). The 
primary limit (Eq.(3.9)) has only powerdaw solutions at infinity unless (3 — 0. We thus tentatively choose 
this value (anticipating later corrections for large but finite N). In this case the solution to Eq.(3.9) is given 
implicitly by 

Uoo - logUoo = y + 1, (3.10) 
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using the boundary condition u(0) = 1. For small y, this equation yields two possible behaviors for u(y), 



«(»)«1±^- (3.H) 

We choose the minus solution to obtain a solution decaying as exp(— y) for large y. The value /3 = 
corresponds to ( = 0, so that the interface remains flat to 0(e) at N = oo. A more careful analysis yields 
the detailed form of the corrections. 

To expand towards finite N, the natural first step is to examine perturbatively the effect of the terms 
dropped in Eq.(3.9). Since the zeroth order solution is exponential for large y, it is immediately clear that 
perturbation theory breaks down in this regime, due to the presence of the term 

Hyu" S> u for y 1//J, 



with u ~ exp(— y). 



The solution Eq.(3.10) is therefore valid (even approximately) only for y <C l//x, and we thus have a boundary 
layer for y l//i. 

Fortunately, we can proceed with the analysis by noting that for y 1, the non-linear terms in Eq.(3.6) 
become negligible, and the equation can be reduced to the form 

u + /3yu' + v! + iiyu" = 0, (3.13) 

where we have allowed for j3 ^ (but 0([i)), which will turn out to be the case for for N large. Eq.(3.13) is 
valid within the boundary layer. It is important to note that both this equation for the tail and the primary 
equation (Eq.(3.9)) are valid in the (asymptotically infinite) region 1 <C y <C 1/fi, which makes it possible 
to match the solutions of the two equations in this domain (see Fig. [2]). 

The linear tail equation (Eq.(3.13)) is second order, and therefore has two independent solutions. For 
Hy,Py <C 1, the solutions are 

u a (y)~ Ce~ v , 

, (3-14) 
u b (y) ~ Cy- 1 ^. 

From the behavior of the primary solution, Eq.(3.10), in the matching region j/ > 1, we see that we must 
choose the solution u a (y) (with possibly a small admixture of Ub vanishing in the N — > oo limit). For a 
given value of (3, both the primary and tail solutions are thus completely determined, giving a uniformly 
valid solution to the full equation (Eq.(3.6)) for large N. 

For short-range correlated disorder, the bare unrenormalized function u(y) decays exponentially (or 
more rapidly) for large y. It is straightforward to see that this exponential decay is preserved by the flows, 
from Eq.(3.6), or directly from Eq.(3.1). Higher order terms from higher-loops can generate at most power- 
law corrections to the initial exponential behavior (from terms with all but one of the R s evaluated at 
= 0). The non-renormalization of the exponential behavior can be seen in a schematic way directly from 
the Hamiltonian. When a particular fast mode is integrated out, as in the previous section, the integral 
performed is of the form 

f [d0 > ]cxp[-0 2 > /2 + .R(0 < +</>>)] . (3.15) 

For large 4> < , the contributions from small and large 4> y can be easily estimated. For small <p > , the argument 
of R is large, so that it can be approximated by an exponential decay R ~ exp(— 0^/2), and the integral 
yields an exponentially decaying function of <p < . For large </> > , the integral is dominated by > w —</><, so 
that now the quadratic term dominates and again yields an exponential function of cj> < , of the exact form 
that corresponds to u(y) ~ exp(— y). 

To be a valid fixed-point function for short-range correlated disorder, therefore, u(y) must have an 
exponential tail at large y. By analyzing the tail equation (Eq.(3.13)) for y ^> 1/(3, we again find two 
possible behaviors: 

u c (y) ~ Cy^-^e-W^, 
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For most values of /3, the required solution u a (y) in the intermediate region 1 <C y <C 1//J, 1//3 will be a linear 
combination of u c (y) and Ud(y), and thus will have a power-law tail at large y. For some special values of 
(3, however, u a (y) will correspond exactly to u c (y), so that the power-law tail of Ud(y) does not contribute. 
This will be the eigenvalue-like condition that determines (3, and hence the short-range roughness exponent 

(su- 
it is, in fact, simple to guess one such value of (3. From the fact that (3 vanishes for N = oo (/_t = 0), 
one expects that (3 = 0(/j,). It is easy to check that for (3 = /j, an exact solution of Eq.(3.13) satisfying the 
matching conditions is u = e~ x ~ y . This value of (3 corresponds to a roughening exponent of 

C«e/(4 + A0, (3.17) 

which we anticipate will be valid for large N. [21] Note that this agrees with a recent replica symmetry breaking 
calculation, which was claimed to be valid in the large N limit. [1] 

To compute the next order corrections to the result (Eq.(3.17)), we need to consider in detail the effects 
of the neglected terms. For the primary solution (y <C 1/m)j these can be computed perturbatively in \i. 
However, they will not affect [3 unless the non-linearities for i/ » 1 are also taken into account, since the 
solution found above can be scaled to match the corrections to the primary solution. We therefore first 
consider the effects of non-linearities for y >• 1. 

To look for solutions close to the original one (m^, in Eq.(3.10)), we make the change of variables 



= exp (^-y - J dy'a(y')j , 



where the final change of variables was made to concentrate on the change of character of the solution for 
y <~ l//_i. With some algebra, Eq.(3.6) is transformed to 

l^ja-b + a 2 - ^+g( T1 ) = 0, (3.19) 

where 

We anticipate that a(t]) will be exponentially small (in for rj ~ 1. Then we can ignore terms of 
order a 2 and also terms of order <j(ri) exp(— rj/ii). We see, however, that g(rf) has a part which depends on 
the behavior of 17(77) for 77 of order l//j (i.e. y of order one) via u, Eq.(3.18). We can obtain this by perturbing 
in (3 and \x about the solution which is good in this region. Anticipating that b <C 1, a perturbative 
calculation performed in appendix A yields: 

w^Uoo-f/iTii ^ e- 1 ^ (I - [1 + Oifj 2 )) fort/»l. (3.21) 

This expression is valid for 77 -c 1 but we will see that the corrections embodied in Eq.(3.18) will be small 
for the desired full solution. From Eq.(3.21), the g term becomes 



e -l-r)/p 

77 



2n 2 

-i- + 1 - 47/ - 2t7 2 + O(n) 
I- 1 



(3.22) 



Neglecting the a 2 (and crexp(— t]/fj,)) terms in Eq.(3.19), we obtain a linear homogeneous equation for 
a. This can be solved straightforwardly by introducing the integrating factor 



F(t7) = exp 



■i/v*- J 



exp \ ![(&- 1)77 + log 77] \ . (3. 28) 
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whence 



a iv) " { c+ Jo df i m [g[f]) b] } • (3 - 24) 

For 77 — > 0, this must match onto the solution Eq.(3.21) implying that, since F(rj) vanishes for small r], the 
integration constant C must be zero. At the other end, since F(rj) also vanishes for large i], a will diverge 
unless the integral in Eq.(3.24) is zero. We thus obtain the integral condition 

/>OC 

/ dr] [b- g^e-^ ^ = 0. (3.25) 
Jo 

With the condition of Eq.(3.25), Eq.(3.24) yields 



dnm ~—b + w^w + 0{ri ) ' (3 - 26) 



for r\ >• 1 and hence 



u(y) ~ y ^(i- 6 ) 2 exp 



-y i 



(3.27) 



for y 3> Up to corrections of order b 2 in cr (which arise from the neglected terms in Eq.(3.19)) this 

agrees with the behavior of u c (y) from Eq.(3.16). We have thus found the desired exponentially decaying 
full solution valid for the full range of y.[22] 

As jU — > 0, both terms in Eq.(3.25) can be evaluated by steepest descents. Expansion of each term 
around the saddle points rj* g {1, 1/2} yields 

1 1 1 

b= — 2"m-[1 -2h + 0{h 2 )], (3.28) 

which is exponentially small for small /j,, justifying our approximations. We are now in a position to obtain 
the roughness exponent from Eq.(3.28). For large N, we have 



j_(i + i 2 -mW 

AT + 4 1 4e + 4 



iV + 2 



(3.29) 

+ ... 



Note that an approximate analysis of the RG flow (Eq.(3.1)) by Natterman and Leschhorn[19] also gave the 
same prefactor, but a different exponentially small correction. If we naively truncate the series after the "1" 
in the last bracket, we obtain for N = 1 

({N = 1) = | (1 + 0.0585) w 0.2117e, (3.30) 

It is clear from the increasing nature of the next term that the series is asymptotic; nevertheless, the 
magnitude of the correction to e/5 compares fairly well with the direct numerical solution of Eq.(3.1) for 
N = 1, which yielded C = 0.2083e. 

In appendix D, we show that in addition to the fixed point found in this section, there is a discrete 
infinite series of fixed points with u(y) changing sign but still decaying rapidly. At this point, whether these 
are physically meaningful is unclear. 
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IV. Stability Analysis 

In this section, we analyze the stability of the short-range fixed point found in the previous section to 
perturbations with both short and long-range correlations. Since this fixed point represents a phase of the 
system, rather than a critical point, there should be no short-range perturbations which are relevant. From 
the scale invariance of the fixed-point equation, it is clear, however, that there is at least a marginal operator 
connecting the line of short-range fixed points with different u(0). This does not alter the physics, since it 
just corresponds to redefinitions of dimensional quantities and is thus an uninteresting redundant operator. 
Intuitively, one expects perturbations with sufficiently long-range correlations to be relevant, causing the 
system to flow to an appropriate long-range fixed point; we will see that this is indeed the case. 

We look for the eigenoperators in the usual way, by considering a function u(y) initially very close to a 
fixed point solution u*(y), i.e. 

u(y)=u*(y)+v(y), (4.1) 

where v(y) is a small perturbation. Inserting this into the flow equation, Eq.(3.6), and keeping only terms 
first order in v(y) gives 

=v + v'(l - u*) + u*'(v(0) -v) + Pyv' + w(v" + u*"v(0)) 

+ 11 [6yu*'v' + 2y 2 {u*'v" + u*"v') + y{u*"{v - v{0)) + (u* - l)v")] , 

where we have chosen to perturb around the fixed point with u*(0) = 1. Note that it is not permitted at 
this point to choose v(0), since it may not be a constant. The right hand side of Eq.(4.1) can be thought 
of as a (nonlocal) linear operator acting on v(y). Just as in finite-dimensional RGs, solutions with simple 
exponential t dependence can be found (The question of completeness is discussed in the footnote [23].) if 
the spatially dependent part obeys the eigenvalue equation: 

v + v'(l - u) + u'(l -v) + (iyv + ny(v" + u") + fiO(uv) = Xv, (4.3) 

where we have dropped the asterisk on u(y) and chosen v(0) = 1 : since eigenvectors are defined only up to 
a constant, we have the freedom to choose a scale for v(y). It is straightforward to show that there are no 
eigenfunctions with v(0) = 0, the only choice not equivalent by a choice of scale to v(0) = 1. This is shown 
in appendix B. The terms in the square brackets in Eq.(4.2) are of order u.uv. For small u,, they are small 
for all y, and will be neglected in what follows. 

As for the fixed point equation, the solutions to the eigenvalue equation can be found in two regions 
and the pieces matched asymptotically. For y <C 1//^, 1//3, the equation can be rewritten in the standard 
form for first order linear differential equations by dropping the 0(/j,uv) terms, 

v'(y) + a(y)v(y) = b(y), (4.4) 

where 



1 - A - u'(y) 
1 - u(y) 



a (y) 

1 - u{y) 

Using the fact that u(y) satisfies Eq.(3.9) (with (5 = 0), and taking some care due to the singularity in u'(y) 
for small y, one finds the solution 

1 M 1_A - U „, 

for y <C l//i, 1//3. In the matching regime, y ^> 1, this solution behaves like 

v - e- x - y [e Xy - 1] /A. (4.7) 
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For y » 1, Eq.(4.3) reduces to a linear equation very similar to the fixed point equation (Eq.(3.13)) in 
this regime, 

Hyv" + (1 + 0y)v' + (1 - \)v = -u' - nyu". (4.8) 

As in the previous section (for the behavior as a function of f3), one expects a discrete set of exponential 
solutions which match onto Eq.(4.7). A natural substitution is thus v — e~ v w, which yields 

fj,yw" + (1 - liy)w' - Xw = (1 - liy)e~ x , (4.9) 

using [3 — ii and u = e~ 1 ~ v , which are valid for small fi. Expanding the solution in a power series, 
w (y) = Em w mV m , results in the following recursion relations for the set {w m }: 

e _1 + Xw 0l 
— /ie _1 A + fi 

2(TT^) + 2(i + 7) Wu ( 4 - 10 ) 

A + [irn 
(to + 1)(1 + /im) 

If the series does not terminate, then the recursion relation at high order simplifies to w m+ \ w w m /m, 
so that w m <~ 1/to!. This implies that for large y, w(y) ~ e y (times a power law of y arising from the 
corrections to the recurrence relation), so that the corresponding cigcnfunction v(y) has power-law decay. 
This implies that the desired short-range eigenfunctions form a discrete set, corresponding to the condition 
that the series terminate at the m th order, for m = 1, 2, . . .. 

For to > 2, these conditions yield the eigenvalues A = — n[i. For to = 1, the matching conditions fix 
wo- We guess that |A| <C 1, and check for self consistency. Under this condition Eq.(4.7) becomes w ~ e _1 y, 
so that Wo = and W\ — e -1 , which satisfies the first recursion relation, as it must for the solutions to 
match, and yields the condition A = in order for w 2 to vanish. This is thus just the redundant marginal 
operator resulting from the choice of normalization of the fixed point solution mentioned earlier. The full 
set of physical short-range eigenvalues is therefore 

A = -2/i,-3M, -Afi... (4.11) 

for small fj,. Note that there is no cigcnfunction with A = — \i. It is absent because of the inhomogeneous 
terms in Eq.(4.9), as discussed above. 



W\ = 

w 2 = 
w m+ i = 



V. Long-Range Correlated Disorder 

We now analyze the behavior for random potentials with long-range power law correlations in cj>. We 
first consider the stability of the short-range fixed point analyzed above to such long-range correlated per- 
turbations. 

The behavior of the power-law eigenfunctions is much simpler than the exponentially decaying solutions 
found above. For any value of A not in the discrete set of Eq.(4.11), the solution which is well behaved at the 
origin has power-law decay at infinity. There is thus a continuum of power-law eigenfunctions with v <~ y~ T , 
corresponding to R(<p) ~ 4>~ J , with T = 1 + 7/2. The associated eigenvalues are 

A = l-/3r, (5.1) 

for all T except those corresponding to the short-range solutions (Eq.(4.11)), i.e. for 

1 (1-2/x) (l-3p) (1-4/x) 
1 ^ 0' ' • ' { ' 

From Eq.(5.1), we see that perturbations are irrelevant for T > l//3.[23] 

We now show that the behavior for distributions with long-range correlations is, in fact, rather simple 
and general. As seen in section III, most solutions of the fixed point equation have power-law tails. It is only 
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for certain special values of £, corresponding to the short-range eigenvalues, that the well-behaved solution 
at y — connects to an exponentially decaying solution at infinity. In all other cases, the solution has a 
power-law decay for large y, which is dictated entirely by the value of £ and e. 

A simple large-?/ analysis of Eq.(3.6) gives the value of ( quoted in the introduction 

Cifl = c/(4 + 7), (5-3) 

for 

(V(<t>, x)V(0', x'))c ~ |0 - 0T^ d (x - x'). (5.4) 

It is easy to sec that this result will be unaffected by the higher-order terms in an RG expansion, since 
higher powers of R(4>) always appear with two derivatives, and multiplying negative power-laws results in a 
more negative power- law. (We restrict our attention to 7 > —2, which is needed to make the problem well 
defined.) Thus, at least within the perturbative RG, only the scale change terms in the RG flows (i.e. those 
in Eq.(3.1) that arc multiplied by e and () arc needed to fix (. The exponent 7 thus fixes the roughness 
exponent at the long-range fixed point exactly. We believe that this result should be strictly true in all 
dimensions, but an actual non-perturbative proof would clearly be desirable. 

For Eq.(5.1), we see that the condition that the long-range correlations dominate and that Eq.(5.3) 
applies is that 7 < 1/Psr, implying 7 < 7 C = e/(sR~ 4. Since this also arises from just the rescaling part of 
the RG flows, we expect it to be true in all dimensions less than four. As 7 decreases, the short-range fixed 
point will become unstable when 7 = -f c ; and for 7 < j c , ( will vary continuously away from (sr according 
to Eq.(5.3). 

A special case of long-range correlated randomness corresponds to interfaces in random field systems 
which have N = 1 and 7 = —1. The general result of Eq.(5.3) implies that (rf = (4 — d)/3 as obtained by 
many authors. [6-8] 

Since the short-range fixed point becomes linearly unstable for 7 below j c , it should be possible to 
observe the instability of the long-range fixed point in the opposite regime, as 7 increases to j c . This is 
actually a subtle problem. If one proceeds with a naive calculation of the cigenfunctions around the long- 
range fixed point, one arrives again at Eq.(4.3). In this case, however, a simple argument demonstrates that 
all the solutions have power-law form. For large y, Eq.(4.3) becomes 

fiyv" + (1 + py)v' + (1 - X)v « AT [1 - M (r + 1)] y- T ~\ (5.5) 

using u*(y) ~ Ay~ r for large y. Since the power-law dictated by the fixed-point function appears on 
the right-hand side, the homogeneous terms must balance this, and the only way they can do so is to 
develop power-law tails themselves. Therefore, all the eigenfunctions have power-law decay at large y. This 
means that the growth of short-range correlations does not manifest itself in the usual manner as a relevant 
eigenvalue. Although we have not calculated this explicitly within the framework of the large N RG, wc 
expect the following behavior: as 7 increases towards j c , the fixed point correlation function of the random 
potential will become more and more like the short-range fixed point with the regime in which the power 
law tail appears moving out to larger and larger y, eventually disappearing for 7 > 7 C .[24] 
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VI. Conclusions 



In this paper we have studied the problem of an oriented manifold with d internal and N transverse 
dimensions. For d near 4, it became possible to treat the zero temperature fixed point in an expansion in 
e = 4— d, in terms of a second-order non- linear differential equation. At N = oo, this equation could be solved 
exactly for the most interesting case of short-range disorder, yielding a roughness exponent (sr = e/(4 + A). 
For large but finite N, boundary-layer techniques were employed to estimate the leading corrections, which 
were found to be non-analytic, vanishing as 2~ N . The magnitude of this non-analytic correction for N = 1 
is comparable to the correction to e/5 found numerically: ( = 0.2083e. The behavior in large N near 
four dimensions saturates a lower bound ( > (4 — d)/ (4 + A) that is believed to be exact. [25] The leading 
correction to this that we have found is indeed positive as it should be. Formally, preliminary analysis of 
our RG flow equations in the opposite limit yields ((A = 0) = e/4. This is believed to be an exact upper 
bound for general A [25]. All known numerical and analytical results do indeed lie in the range 

4 — d A A-d ,„.,. 
<C<^- (6.1) 



4 + A 



It is interesting that our large N e-expansion result agrees with the large A results of Mezard and 
Parisi, claimed to be valid for general dimension (2 < d < 4) by replica symmetry breaking techniques. 
A preliminary large A analysis of the zero temperature minimization problem corresponding to our RG 
procedure does not suggest that the 0(e) result should be exact for e < 2 in this limit, although further 
study may change this conclusion. The analysis does, however, suggest that the next corrections to (, which 
are probably of 0(e 3/2 ), may be calculable without taking into account the effects of multiple local minima. 

Both extensions of the e-expansion and investigation of the large A limit beyond the e-expansion are 
worthwhile future endeavors. An important question is whether or not the "replica-symmetry breaking" used 
in the variational ansatz by Mezard and Parisi has any well-defined physical interpretation beyond that of 
the general scaling picture of manifolds in random media (discussed in detail in the directed polymer context 
by Fisher and Huse[4]). Answering this might bear fruit for understanding other random systems, such as 
spin and vortex glasses. Finally, it is possible that some of the techniques used here might be applicable for 
other problems such as periodic (e.g. charge density wave or flux lattices) or non-periodic (e.g. polymerized 
membranes) elastic manifolds which exist in random media of the same dimension. 
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Appendix A. Perturbative Calculation in Primary Region 

For the primary region, y <C simple perturbation theory can be used to find the effects of the O(u) 
terms of Eq.(3.6). The solution is expanded in a power series in ^, 

u = Uoc + A*tii + H 2 u 2 + . . . , (A.l) 

and terms are grouped order by order, assuming (as confirmed in section III) that (3 = O(u). The zeroth 
order equation is just Eq.(3.9), while the first order terms give 

- «oo) + Ul(l - O = iV^'oof + 2j/Voc<c + J/Ooo - 1) - (P/riVv'oo, (A.2) 

where u oa (y) is the solution of Eq.(3.9). The right hand side of this equation can be rewritten completely in 
terms of by using Eq.(3.10) to eliminate the y dependence, and Eq.(3.9) to eliminate derivative terms. 
After some lengthy algebra, one finds 

ui(l-« 00 )+«i(l-«^) = fr(u 00 ) ) (A.3) 

where 

H(x) = H (x) + H (x), (A.4) 
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with 

H (x + 1) = [-2/x 2 - 2/x + 3x + 3] + log(a; + 1) [4/x 3 + 6/x 2 -l/x-3] 
+ log 2 (.x + 1) [-2/x 4 - 4/x 3 - 2/x 2 ] , 

and Hp{x) = [x - log(x) - log(x)/(x — 1)] . 

M 

Since Eq.(A.3) now has no explicit y dependence, we switch to the dependent variable Moo, using 

dUl !toc 



ay duoo 



du n 



(A.5) 



(A.6) 



Eq.(A.3) becomes 



which has the solution 



du\ 

dUoc 



+ 



Ul 



«oo(«oo - 1) 



Ml(tioo) 



Uqo f 
1 - Uoo Ju 



-H(Uoo)/u 

1 1-x 

— T7—H(x)dx, 
x A 



(A.7) 



(A.8) 



where the boundary condition u\{y = 0) = has been imposed. The important limit for the matching 
carried out in section 3 is y » 1, corresponding to ~ exp(— 1 — y) — > 0. In this limit, the integral yields 



IMixiy) w e 



(A.9) 



In section III, it was found that |/3 — <SC /U, so that the first term in the brackets can be neglected. The 
full solution is then still a pure exponential in the matching region, but with a different coefficient: 



i(i/)^" H (l-/i) forl«y«l//i. 



(A.10) 



Appendix B. Eigenfunctions with v(0) = 

By choosing v(0) = 0, we arrive at an equation for this case analogous to Eq.(4.3), 

v + (1 - u)v' - it't> + 0yv' + Liy(v" + u") = Xv, (B.l) 

where we have already neglected the [iO{uv) terms. In the perturbative region (y <ti the equation 

analogous to Eq.(4.4) is 

v'(l-u) + (l-\-u')v = 0, (B.2) 

which in this case is homogeneous, and correspondingly simpler to solve. Using Eq.(3.9)and Eq.(3.10), one 
finds the general solution 

v = Cu 1 - x /(l -u). (B.3) 

As y — > 0, u — > 1 — \/2y, so that 

v(y) - -g=. (B.4) 



For u to be non-singular for small y, the only choice is C — 0, so that no such eigenfunctions exist. 

This is in some respects a surprising result, since for any perturbation of u*(y), one should be able to 
rescale the resulting function to leave u(0) unchanged. Therefore one might expect the stability analysis 
to be phrased precisely in terms of those perturbations which do not change u(0). The rescaling, however, 
is equivalent to adding some amount of the marginal eigenfunction which moves along the fixed line. The 
function representing the combined effects of perturbation and rescaling (with v(Q) = 0), is, in this case, not 
an eigenfunction. 
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Appendix C. Schematic Minimization 

To understand the effects of multiple minima in the random potential on the validity of the RG used here, 
it is instructive to consider a simple model problem representing the iterative minimization at a particular 
length scale. This provides a physical derivation of the renormalization group flow equation (Eq.(2.21)). 
Schematically, we imagine integrating out a single Fourier mode of momentum p. At zero temperature, 
this reduces to the problem of minimization over this TV— dimensional vector. The renormalized potential 
will then be given by 

Vfl(^ < )=min|^ 2 > +V(0 < +0 > )|, (C.l) 

where the magnitude of the momentum cut-off p = A has been set to one for convenience, and, more 
importantly, we have ignored the spatial dependence of the potential V(0, x). A similar formulation of the 
RG would appear in the treatment of elastic manifolds on a hierarchical lattice. Although Eq.(C.l) may 
appear to be an unreasonable approximation, we will find that the lowest RG order flow equations that we 
have used in this paper are exactly reproduced. Indeed, the treatment of a suitably modified version of 
Eq.(C.l) which takes into account the spatial dependence does not differ substantially from the approximate 
version described here (though the required notation makes it rather more cumbersome). 

For the purposes of simplicity, we will concentrate on the case of N = 1, which, for small e, does not 
differ much from the case of general N. For the remaining part of this appendix, we introduce the notation 
x = 4> < and y = <py. The N = 1 problem is then 

e 1/2 V R (x) = min U(x, y) = min ( \y 2 + e^V ( x + y )\ ( C .2) 
v v \ 2 J 

where a factor of e 1 / 2 has been extracted to make V of order 1. To estimate quantities, we use the fixed-point 
values for the correlation function of V(y), 

(V(x)V(x')} =R(x - x 1 ) = R(x - x')/e, 
(V(x)) =0. 

Because the minima in Eq.(C.2) will be at small y for small e, it is the small-distance behavior of R(x) which 
will determine the behavior of the model. Based on the behavior of the fixed point function R(x) in the RG 
of section II, we assume that R(x) can have a discontinuity in its third derivative at x = 0. 
The extremal condition for Eq.(C.2) is 

y = e 1 ' 2 F(x + y), (C.4) 

where 

F(y) ee -V\y) (C.5) 

is the force at the "position" y, and primes have been introduced to denote derivatives. From Eq.(C.3) the 
correlations of the force at short distances are 

(F (xf) =1, 
((F(x)-F(0)f)~\x\, 

for \x\ < 1. 

The perturbative RG performed in section II is equivalent to assuming a perturbation series for y in e 1 / 2 
in Eq.(C.4) and expanding F(x + y) to obtain solutions order by order in e 1 / 2 . There are two possible ways 
in which such an expansion may break down. Most obviously, the linear behavior in Eq.(C.6) implies that 
derivatives of F(y) are typically infinite, so that the analyticity assumed by a power series in e 1//2 may break 
down. Secondly, the extremal condition (Eq.(C.4)) may have multiple solutions, and an iterative solution 
may converge to any of these, including both maxima and non-global minima. 

Note that if R(x) were smooth, then ([V"(x)] 2 ) < oo and the full potential U in Eq.(C.2) would have 
strictly positive curvature with high probability for small e and hence a unique minimum. The apparent 
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(C.3) 



(C.6) 



assumption of analyticity of R(x) may be circumvented by the use of a different iterative procedure to 
find the minimum of the potential U . The simplest such method, which reproduces the results of the 
perturbation series for the analytic case, is a version of gradient descent (see Fig. [3]) One defines a sequence 
of approximants {yo, yi, yi, . . .} to Eq.(C.4) by 

y n+1 = e 1/2 F(x + 2M), (C.7) 



for n = 0,1,2,.. . 

It is interesting to note that the stability properties of such a mapping discriminate between minima 
and maxima. Letting y n = y* + Sy and linearizing, one finds 

Sy n+ i = -V"(y*)Sy n , (C.8) 

which is stable for |V"(y*)| < 1. Since the curvature of the potential U in Eq.(C.2) is just 1 + V"(y), 
this condition excludes all maxima (as well as minima which are sufficiently narrow). Note, however, that 
the iterative scheme does not guarantee convergence, and one cannot rule out limit cycles or other more 
complicated behavior for particular realizations of the disorder. Furthermore, even if it docs converge, it 
may not be to the desired global minimum. As we shall see, however, we will obtain an estimate of the 
global minimum with sufficient accuracy for our present purposes. 
Iterating Eq.(C.7) yields the first few approximants, 

2/0 =0, 

Vl =£1/2F(X) ' (C 9) 

y 2 =eV 2 F(x + eV 2 F(x)), 

etc 

Using Eq.(C.6) the corrections at each level of approximation may be estimated. From Eq.(C.7), 

Vn+i - Vn =e 1/2 [F(x + y n ) - F(x + y n _i)] 

1/2, ,1/2 {^.LU) 

~C ' \Vn - Vn-l\ ' , 

with a random coefficient of 0(1). Iterating this, an infinite series of non-trivial powers of e appear. For the 
n th approximant, we thus have 



y 



n « Ae 1 ' 2 + Be 3 / 4 + Ce 7 / 8 + ... + Ze^-^ 2 ". (C.ll) 



The difference from the perturbation series result for the analytic case appears first in y 2 , via the appearance 
of the e 3 / 4 term in Eq.(C.ll). 

From Eq.(C.4), we may find an upper bound for the size of the region in which minima are likely to 
exist. (In fact, we only calculate this bound for the region within which there are extrema. For N = 1, the 
furthest out extrema are in fact always minima, but for large N this difference may be important.) Suppose 
one extrema is located at the point y*, satisfying Eq.(C.4). Then for another extremum to be located within 
a distance Sy, the condition must again be satisfied at the point y + Sy. For large Sy, this is clearly extremely 
unlikely, since the linear term grows while the random force remains bounded and of 0(e 1 ^ 2 ) with high 
probability. For small Sy, the variations of the force grow like (Sy) 1 / 2 from Eq.(C.6), i.e. faster than the 
linear term in Eq.(C.4). Thus there will be a length scale below which the variations of the force dominate, 
and other extrema are possible. Equating the two terms in Eq.(C.4) gives 

Sy^e^ 2 \Sy\^ 2 ^Sy<0(e), (C.12) 



which is the desired upper bound on the separation of extrema. We see that the separation between extrema 
is smaller than any of the correction terms obtained in (C.ll). This is illustrated in Fig. [3]. This suggests 
that the corrections due to multiple minima appear at higher order in e than the iterative corrections which 
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arise from the non-analyticity of V(y). To check this, we must analyze how the corrections to y affect the 
renormalized potential and its correlations. 

To estimate the corrections from the terms in Eq.(C.ll), it is useful to write the random potential is 
the form 

V(x + y) = V(x) + yV\x) + W(x; y), (C.13) 



where 



ry 

W(x;y) = / [V'(x + z)-V'(x)]dz. 
Jo 



(C.14) 



The renormalized potential can then be evaluated by inserting the iterative solution Eq.(C.ll) into Eq.(C.2), 
yielding 



e^V R (x)=e^V(x) + ^e+^e^ _ ; 



+ C 2 {x) c7/4 + 



(C.15) 



+ A(x)B(x)e 5 ^ + A{x)C{x)e 11 / & + ■■■ + e^ 2 V\x) A(x)e^ 2 + B(x)e 3 ' i + ■ 
+ e 1 / 2 W[x;A{x)e 1 / 2 + B(x)e 3 / i + •••]. 

This expression simplifies somewhat when the function A(x) = F(x) = —V'{x) from Eq.(C.9) is inserted, 



e^ 2 V R (x) =e 1 ' 2 V{x) - ^6 + ^M e 3/2 + . . . 

+ B{x)C{x)e 13/H + ■■■ + e^ 2 W[A{x)e^ 2 + B(x)e 3 / 4 + •••]. 



(C.16) 



The fact that the second term is negative reflects the approach to the minimum. Note that the cross terms 
AC, AD, etc. have now canceled. The RG flow requires renormalization of the correlation function of the 
disorder, Eq.(C.3). The renormalized potential in Eq.(C.16) will have a non-zero expectation value, since the 
minimization procedure decreases the energy for all realizations of V(y). To find the renormalized correlation 
function, therefore, it is necessary to take the truncated (cumulant) expectation value, 



eR R {x)=e(V R {x)V R (0)} c 
=e(V(x)V(0)) 



l-e 3 / 2 (A 2 (x)V(0)+A 2 (0)V(x)) c + l-e 2 (B 2 (x)V(0) + B 2 (0)V(x)) c 



+ -e 2 (A 2 (x)A 2 (0))c + 2e(V(x)W[0]} c e + e(W[x]W[0]}c 
-le 3 l 2 {A 2 {x)W[Q]+A 2 {Q)W[x])c + \e 2 {B 2 {x)W[Q]+B 2 



(C.17) 



where in this expression W[x\ means the full expression from the last term of Eq.(C.16), evaluating all 
internal coefficients (i.e. A(x), B{x) etc.) at the point x, and we have dropped all terms of explicitly higher 
order than e 2 . The function W[x] is, however, itself small, so that some further terms can be dropped. From 
Using Eq.(C14) and the fact that F{x) - F(x + z) ~ sfz (Eq.(C.6)), we see that W(x;y) ~ |y| 3 / 2 . Since 
y <~ e 1 / 2 , W[x] ~ e 3 / 4 + O(e), where the O(e) term arises from the B(x) in Eq.(C.ll) and higher order terms 
in A(x); this O(e) term is needed to obtain the V R correlations to 0(e 2 ) but higher order terms are not. In 
addition we see that of all the terms involving W in Eq.(C.17), only the VW term will contribute at this 
order. We thus see that to 0(e 2 ), we may drop all terms in y beyond the B(x) term. To this order, we may 
thus use the second order iterative solution 



e^VRix) w e^V \x + e x l 2 F[x + e^ 2 F(x) 



1 

+ r 



F[x- 



,1/2 



F(x) 



(CAS) 



The correlations of V R (x) and V R (0) can be calculated directly from this form of V R yielding 



Rr(x) - e(V R (x)V R (0)) c = R(x) + - [R"(x)f - R"(x)R"(x) + 0(e 5 / 2 ). 



(C.19) 
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To this order, the correct answer can be obtained from Eq.(C.17) by expanding W(x: y) formally as W(x; y) w 

2 

HfV"(x) and averaging directly the WV term in Eq.(C.17). This is valid because only one V" appears here. 
To analyze the WW term, an expansion in y fails and (WW) <~ e 3 / 2 rather than the naive e 2 . We thus expect 
that the effects of the nonanalyticity of R will affect Rr at order e 5 / 2 . These terms need to be balanced by 
adjustments to £, suggesting 0(e 3 / 2 ) corrections to our 0(e) result for (. 

So far, the effects of multiple extrema have not been included. Their effects can be estimated by 
including a further correction term in y, 

y = Ae 1/2 + Be 3/4 + • • • + r]e + ■ ■ ■ (C.20) 

By repeating arguments along the lines of those above, it is a simple matter to estimate the leading corrections 
due to a non-zero 77. One finds that the first contribution to R(x) occurs at 0(e 3 ). This is higher order than 
all the leading corrections from the non-analyticities arising in the iterative procedure. 

Thus, although the singularity in R"'(x) at x — is associated with the existence of many extrema 
(since V"(0) has infinite variance), the direct effects of these multiple extrema only show up at higher order 
in e; to the order needed here, choosing any of the minima provides enough accuracy. 

The simple approximation to the renormalization group flows analyzed in this Appendix suggests that 
there will be corrections to £, starting at 0(e 3 / 2 ), with an apparently infinite sequence of higher order 
corrections appearing before 0(e 2 ), at which order the effects of multiple minima begin to appear. Although 
the picture is quite appealing, the results should not be taken as definitive predictions of the powers involved, 
since a complete analysis should involve a self-consistency condition to determine the small x behavior of 
the correlations. Such an analysis may well involve a boundary layer for small x with smoothing of the 
fixed-point function on scales smaller than some (> 1) power of e. 

It is straightforward to extend the analysis of this Appendix to general (fixed) N in the limit of small e. 
Since in the limit of large N, £ is formally small, even if e is not small, one might hope to be able to justify 
truncation of the RG flows for all e (or at least e < 2) for N large. We have not been able to do this, and, 
indeed, preliminary indications suggest the opposite conclusion: that even for large N higher order terms in 
e are needed. A more detailed study of this limit would clearly be instructive. 



Appendix D. Multicritical Short-Range Fixed Points 

As remarked in section III, Eq.(3.13) possesses a discrete family of solutions which are well-behaved 
at the origin and decay exponentially at infinity. These can be matched onto the primary solution to yield 
additional fixed points of the RG flows. To find these, we perform a power series expansion similar to the 
one used for the stability analysis (Eq.(4.10)). Defining 



u = e ^'w, 



w 



w m y 



(D.l) 



Eq.(3.13) yields a simple recursion relation for the set {w m }: 

'/ 'fi- 1 + fim 



w m+ i 



_(m + l)(l + /zm)_ 

If the series does not terminate, the large m behavior of the coefficients is 

1 



(D.2) 



(D.3) 



so that u(y) decays more slowly than an exponential. A short-ranged u(y) is obtained whenever the series 
terminates, which yields the condition 



0-- 



ft 



1 + mjj, 



0,1,2,.... 



(D.4) 
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In terms of the roughening exponent, 



< = I+WT»Z (D ' 5) 



The case m = corresponds to the simple exponential found in section III, while for higher m the solutions 
have some oscillations and correspond to smaller values of £. 

It is a simple matter to extend the results of section IV to calculate the stability around the new fixed 
points. One finds that for the m th fixed point, there are m relevant eigenvalues corresponding to short-range 
correlated perturbations, so that the solution found in section III is stable, while the remaining solutions 
represent a hierarchy of multicritical solutions. 

Although such solutions exist formally, we have not fully investigated the criteria under which these 
solutions represent physically meaningful fixed points. At least initially, the function R(4>) is highly con- 
strained by the positivity condition for the probability distribution of V(4>). In particular, taking the Fourier 
transform of the potential-potential correlation function, we must have 

V(k)V(-k)^ = R(k) > 0. (D.6) 

Although the interpretation as a correlation function suggests that this positivity property is preserved by 
the RG, the non-locality (in k) of the terms in the RG flows generated by fluctuations has prevented us from 
finding a simple proof. Nevertheless, it seems likely that the positivity is preserved. It is straightforward to 
check the multicritical fixed points obtained above for this criterion. If they do not satisfy Eq.(D.6) they 
cannot be physical. A simple computation for the first multicritical solution, 

^(0)=(l-M 2 -^ex P f- 7 ^V (D.7) 

yields the Fourier transform 



2 ; — v 2(1 + M ) 



= Q + \ (1 + exp (-^-(1 + M )) , (D.8) 

in the large N limit, which satisfies the positivity criterion (Eq.(D.6)) and thus might be physically attainable. 

If we restrict consideration to distributions in which the function u(y) has no zeros (or equivalently 
R(\4>\) has no non-trivial extrema), it is possible to show, however, that the multi-critical solutions, which 
have at least one zero are inaccessible. A simple argument proceeds as follows (see Fig. [4]): Consider the 
evolution of an initial function uo(y) which has no zero crossings. For the function uo(y) to evolve into one of 
these multicritical solutions, it must at some intermediate stage when it first has a zero be tangent with the 
y axis at some point t/j (like the function Ui{y) in Fig. [4]). (Note that since the behavior at 00 is preserved 
by the flows, the zero cannot come in from 00 and avoid the tangency condition.) At this intermediate point, 
the function must obey 

Ui(yi;ti) =0, 

«J(w;*i)=o, (D-9) 

and u"(yi,ti) >0. 



From Eq.(3.6), one then finds 

du(yf,t) 



dt 



Huo(t i )y i u"(y i ,t i )>0, (D.10) 



so that the putative point disappears. Therefore such a zero cannot occur, and any function u(y) with zeros 
is inaccessible from an initial u without zeros. 

At this point it is unclear whether the formal multicritical fixed points found here are accessible for 
less restrictive initial correlation functions, and, if so, what is their physical significance. In particular, one 
might expect the m = 1 critical point to separate two phases with different behavior. If one of these is the 
rough phase analyzed in this paper, what is the nature of the other phase? We leave these as intriguing open 
questions. 
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Figure Captions 



Fig. 1: Diagrams representing the generation and feedback of the three-replica term S (see the third line of 
Eq.(2.7)). Fig. [la] shows how such an operator is generated at second order in R by terms with one 
contraction. Since the momentum of the internal line is within the momentum shell, the diagram 
only contributes at large momentum. Summing up all such terms resulting from the expansion of R 
(Eq.(2.22)) gives Eq.(2.18). Since the fixed-point value of S is 0(e 2 ), the only potentially dangerous 
contribution for our analysis is first order in S. This feedback comes from diagrams such as that of 
Fig. [lb], with a single loop. Since the original S vertices were generated only at high momentum, 
such terms do not renormalize R. 

Fig. 2: Regions of validity of the primary (perturbative) solution, Eq.(3.9), and the tail (linearized) solution, 
Eq.(3.13). In the large N (small /i) limit, the size of the matching region (1 <C y <SC l//u) grows 
without bound. 

Fig. 3: Graphical illustration of the iterative minimization of Eq.(C.7). Given a guess yi for the location 
of the minima, the next approximation is found by following a vertical line at this value of y until 
it intersects the random force curve. Extending a horizontal line to the 45° line through the origin 
(representing the uniform restoring force of the harmonic potential) gives the value of y for the next 
iteration. The second two approximants, y\ and y 2 , resulting from the initial y = are shown here. 
The many intersections between the random force curve and the 45° line in the figure represent 
multiple extrema which occur on smaller scales (see Eq.(C.20)). 

Fig. 4: Illustration of the preservation of the lack of zeros by the RG flows (Eq.(3.6)) for the function u(y). 

For the initial function uo(y) to develop into the final function Uf(y), with internal zeros, it must 
pass through an intermediate stage Ui(y), at which it is tangent at some point yi with the y axis. 
Eq.(D.lO) shows that such a point of tangency is repelled, so that the putative crossing does not 
occur. 
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